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THE GEMPAK BARNES OBJECTIVE ANALYSIS SCHEME 
S, E. Koch, M. des Jardins, and P. J. Kocin 

ABSTRACT 

GEMPAK is an interactive computer software system developed at the Severe Storms Branch 
of NASA’s Goddard Laboratory for Atmospheric Sciences for the purpose of assimilating, analyzing, 
and displaying various conventional and satellite meteorological data types. The Barnes (1973) ob- 
jective map analysis scheme possesses certain characteristics that allowed it to be adapted to meet 
tlie analysis needs of GEMPAK. Th(ise characteristics and the specific adaptation of the scheme to 
GEMPAK are described hr re. A step-by-step guide for using the GEMPAK Barnes scheme on an 
interactive computer (in “real-time”) to analyze various types of meteorological datasets is also 
presented. 

The GEMPAK Barnes scheme is unique in the way in which it achieves a balance of objectivity, 
versatility, and practicality. Demonstration of th.ese qualities is accomplished by applying the 
scheme to both a high quality, uniformly distributed radiosonde data set, and to a non-ujiiformly 
distributed data set of undetermined quality composed of satellite-derived cloud motion vectors. 
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THE GEMPAK BARNES OBJECTIVE ANALVSIS SCHEME 


L AN OBJECTIVE ANALYSIS SCHEME FOR GEMPAK 

One of the basic program objectives of the Severe Storms Branch (Code 914) of the Goddard 
Laboratory for Atmospheric Sciences (GLAS) is “to study the subsynoptic to mesoscale processes 
which create the favorable environment for the development and evolution of severe convective 
storms by using the recently acquired SESAME data sets, TIROS-N (and VAS) temperature and 
moisture soundings and satellite derived cloud tracked winds” (p. 44 of Simpson et al„ 1980), In 
order to assess the impact that merging of the various data sources has upon the dynamic computa- 
tion and subsequent interpretation of storm-environment processes, the GEMPAK (GEneral 
Meteorological data assimilation, analysis, and display software PACKage) was created. This system 
is designed to (1) convert tlie various types of data to a standard file structure, (2) vertically and 
horizontally interpolate the standardized data to uniformly spaced grid points, (3) perform data 
assimilations to initialize numerical models, (4) compute diagnostic quantities, and (5) display the 
products onto various maps. The GEMPAK has been developed on the AOIPS (Atmospheric and 
Oceanographic Information Processing System) at NASA/Goddard. The adaptation of the Barnes 
(1973) objective analysis scheme to GEMPAK for the purpose of performing the interpolation func- 
tion (step 2) is discussed in this report. 

The term “objective analysis” refers to a process by which data observed at irregularly spaced 
points in space (and/or time) are numerically processed to retrieve tlie two- or three-dimensional 
structure in the spatial distribution of sorrs' -"eorological field parameter. Usually, the observed 
parameter values are computed at a regularly spaced array of grid points in order to provide the 
basis for mapping isopleths of the parameter distribution, and to allow quantitative diagnostic and/or 
predictive calculations to be made. 

Tlie process of objective analysis results in smoothing of both high-frequency, short-wavelength 
features commonly referred to as “noise” and of gross data errors in the data distribution. It is 


necessary to suppress “noise” and reduce error amplitudes significantly, particularly before compu- 
tation of meaningful estimates of such derivative quantities as divergence and vorticity from the 
grid point values, 

The degree of smoothing is governed by the selected v/ciglits applied to the data. The Barnes 
scheme uses a Gaussian weighted-averaging metliod to assign a non-zero weiglit to each observation 
according to its distance S from a mesh point relative to the number and distance of all other obser- 
vations within some “cutoff radius” Rg, beyond which all weiglits are assigned a value of exactly 
zero. Given the Gaussian nature to the scheme, an “influence radius” R can be defined as that S 
where the weight falls to a value of exp (-1) = 0.37 (as discussed in section 2). However, tlie weight 
function is not tied to the chosen value of R^, and in facta// the data points in a given dataset can 
be incorporated in determining each grid point value if desired, by simply setting Rg = «> (norniiilly, 
for reasons of computational economics, Rg is set to a finite value with no loss in analysis accuracy 
as long as Rg is significantly larger than. R, as discussed in section 3e), 

It is of some interest to provide some comparisons between the Barnes (1973) and Cressman 
(1959) objective analysis schemes, as diey botli are weighted-averaging techniques in comm, on usage. 
Perhaps the most important difference is that tlie Cressman weight function is determined by the 
choice of R, since tlie weiglits do not asymptotically approach zero vdth increasing S as they do in 
the Barnes technique, but instead abruptly become zero at S = R(= R^), This nature of the Cress- 
man scheme can presejit serious difficulties when the data distribution is non-uniform. 

All weighted-averaging, “successive correction” (see below) objective analysis methods like 
. those developed by Barnes and Cressman are designed to work best on two-dimensional data fields 
of rather uniform data distribution. Those fields that have very irregular data densiticj over the 
domain can present a problem known as “ballooning,” characterized by large amplitude and phase 
distortions in tlie neigliborhood of any grid point whose value is determined primarily by the value 
of only one inaccurate datum. In general, this problem can be reduced by requiring that at least 
several observations be used in the calculation of each grid point value. In the Barnes technique, 
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this is accomplished by extending tlic distance Rg (with no resulting effect upon the weight function). 
In some applications of the Cressman technique, the current scan R is locally increased to insure that 
a sufficient number of data influence each grid point value (Inman, 1970); however, "noise” intro- 
duced by such a locally varying weiglit function must be suppressed with additional numerical filters, 
thus producing an unknown final response and requiring additional computer time. TJic other cor- 
rectional method used in the Cressman technique is to increase R sufficiently on the first pass through- 
out the entire domain, yet such a procedure means more passes are required to achieve the desired 
final response. 

The Cressman and Banies techniques both employ the method of successive corrections whereby 
an adjustment is made to the first pass analysis by decreasing R in the second pass through the data 
to restore the amplitude of large wavenumber components suppressed in die first interpolation- 
filtering pass. In the Cressman case, neither the number of additional passes nor the value of tiie 
second pass R are governed explicitly by the data distribution, thus neither is the filter response.^ 

An advantage of the Barnes technique is that (in the 1973 version) only two passe.*} are required to 
achieve the desired pattern recognitio n whereas typically (but ratlicr arbitrarily) four to six scans 
are employed with the Cressman technique. 

Classical sampling theory (Peterson and Middleton, 1963) dictates that a wave whose horizontal 
wavelengtii does not exceed at least twice the average observation spacing (2An) cannot be resolved 
since five data points are required to describe a wave and its derivatives. Random errors in the ob- 
servations generate fictitious 2An waves (Barnes, 1964); therefore, it is desirable to filter these from 
the analysis as much as possible. The Barnes (1973) technique has the highly desirable property that 
the weight function constants can be chosen upon the basis of the data distribution prior to the 
analysis so that pattern scales resolvable by the data will be revealed to a known response amplitude. 


Stephens and Stitt (1970) show that an optimum choice for R on the first pass can be made in terms of An 
for a uniform data distribution. However, the optimum R is not well defined theoretically on the second pass, so 
the final filter response is rather arbitrary. 
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To suinmariztj, the Barnes (1973) technique Is selected because of the following advantages 
(listed by him). 

1 . Tile scheme is a computationally simple, Gaussian weiglited-averaging technique which 
assigns a weight solely as a known function of distance between datum and grid point. 

2. Since the weighting function decreases to zero asymptotically, the Infiuence of data may 
be extended any distance to insure that a sufficient number of observations influence each 
grid point value without having any effect upon either the weight function or the response 
clharacteristics. 

3. The low-pass filter response characteristics can be determined prior to the analysis so that 
pattern scales resolvable by the data distribution will be revealed. 

4. Only two passes through the data are required to achieve the desired scale resolution 
because of the rapidity with which convergence is reached (sec appendix). Bven wJien a 
large influence radius (weight foctor k in the next section) k riiosen to reduce noise due to 
variations in observation density, convergence is attainable (because of the numerical 
convergcitce parameter y discussed below). 

5. Small scale “noise*’ is sufficiently filtered from the analysis after only two passes. 

The reader interested only in applying the GEMPAK Barnes sclicme to a particular problem 
can skip the next section without loss of continuUy. 

2. CHARACTERISTICS OF THE BARNES OBJECIIVE ANALYSIS SCHEME 

Barnes 0964) has developed an objective analysis scheme that utilizes a Gaussian wciglit 
function in the spatial domain based on the supposition that two-dimensional data distributions can 
be represented by Fourier integrals. Let represent the distance between the (i,j) grid point and 
an observed datum f(x^ Further, let ic be the weight parameter that determines the shape of 
the filter response function. Then, the weight function w^ is expressed as; 
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W„I * cxp(-r„2/K), (n 

Por a given choice of k, an “innuciico radius” may be thought of as that r,„ wJierc vv,,, « c'** . 
It Is shown later that n can be chosen to reveal those scales resolvable with the particular data 
distribution. 

Barnes (1973) has modified the earlier method in order to decrease the amount of computer 
time necessary to achieve “desired response” at small wavclengllis. This modincation consists of 
applying only a single correction pass through (iteration upon) the initial interpolated lield g„(i,j), 
rather than making several iterations as before. To accomplish tliis, parameter u is decreased from 
its first pass value («^) to its correction pass value of 


«l = T«o« 


by using a “numerical convergence parameter” 7 ( 0 < 7 <!) that forces a high degree of convergence 
(agreement) between the observation field f(x,y) and the correetlon (.second) pa.ss interpolated field 
gj (x.y). Of course, it is not always desirable to have the interpolation field fit the data exactly ; by 
manipulating the value of 7 , the analyst obtains the desired response following the second data pass. 


The first pass yields an interpolated field given by: 


go f(x,y)Do* 


where the spectral response function 


Dq = exp -/^Q(7r/^)2j 


derived by Barnes (1964) is of the form of a low-pass filter that attenuates signals at small horizontal 
wavelength X. Figure 1 shows that use of smaller values of the dimensionless weight parameter 
results in greater filter response during the first pass, particularly for the short waves. 
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Figure 1 . Relationship of first pass response to dimensionless wavelength X^=X/L for various 
values of dimensionless weight parameter w*= kJI?, in which L is an arbitrary scale 
length based upon the data distribution. 
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The “correction (second) pass” grid point field g| (i,j) is the result of adding to the first pass 
field the smoothed residual difference between the observed data values and the first pass estimated 
values “at” the data locations g^fXiy), or 




where (after Barnes, 1973), 


D, 


exp[~Kj (a/X) 

= exp [-7/^0 (a/X)-| 


DT 


is the response function corresponding to the weight function 


Wm == exp(~r„,2/^«j , 


(5) 


i(>) 


(7) 


A simple bilinear interpolation between the values of g^dJ) at tjie four surrounding grid points can 
be used to obtain an estiniute for gj)(x,y) at each ilata location. The aeliial euiTeclion pass value 
“at” each grid point is computed as the sum of the weighted averages from the two passes willi M 
observations according to (after (5)): 


Si 


Ni M 

■£ V w„;tr(x,„,y,„)-So(x,„,y,„)l 
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( 8 ) 


m = l 


Since it follows from (5) that (omitting arguments for convenience) 


D, = 


g] -go 


^ f'-go ’ 


(9j 


the function ( 6 ) miglit be referred to as “a difference field response function.” The true correction 
pass response function (i.e,, that one which corresponds to the first pass response function given by 
(3)) is: 


whicii upon substitution from (5) gives 



D 



D, + (! 


D,)Dj 


( 10 a) 


( 10 b) 


and upon further substitution from ( 6 ) we have (Barnes, 1973) 

of = D„(l+Dj-‘ -Dj), (in 

Function is the proper measure of the degree of analysis convergence, or in ■>ther words, 
how closely the interpolated values agree with the observed ones after a second pass through the 

data. Making a second pass will increase the degree of convergence when 0< y< 1, and particularly 
so wheny< 0,5 (Figure 2). Notice that the greatest increase in response urs at the shorter 

wavelengths (small for a given choice of according to Figure 1). 


A mathematical analysis of the effect of making additional passes (N>2 iterations) through the 
data appears in the appendix. It is proven there that no real benefit can be gained in making more 
than one correction pass because of the rapidity with which convergence is approached when 7 is 
chosen small enough (Figure 2). It is also proven that the 1973 version of the Barnes objective 
analysis technique is absolutely convergent (although the 1964 version has been shown by Barnes 
(1964) to be convpigent also, several more passes are required to reach the same degree of conver- 
gence as with the 1973 version). This fact enables the analyst to control the amount of small-scale 
detail to be revealed in the analyzed data fields. 
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Figure 2. Correction pass response function as a function of first pass response and 
numerical convergence parameter 7 (from Figure 3 of Barnes, 1973). 



These properties of the. Barnes analysis sclieme, (a) weiglit selection based wholly upon the data 
spacing, (b) filtering of the 2An v/ave from the analysis, (c) adequacy of only two data passes to achieve 
desired fit to the data, and (d) insurance that a sufficient number of observations influence each grid 


point value, are all incorporated into the GEMPAK Barnes analysis scheme described in the next 
section, Other factors, such as objective determination of grid size and accounting for variable 
ohservalion ilensity, aie unique rcaturcs of the GEMPAK version of the Barnes scheme. 


3. THE BARNES SCHEME TAILORED TO GEMPAK NEEDS 

Certain factors were considered in accommodating the Barnes Cl?'73) objective analysis sclieme 
to the specific needs of GEMPAK; 

1. An agreeable balance between user manipulability and the objectivity inherent to the 
Barnes scheme must be found, 

2. Differences in the nature of the various conventional and satellite data must be taken Into 
account. In particular, satellite-derived cloud motion wind data tend to occur in clumps 
separated by varying distances because of the problem in finding a sufficient number of 
suitable, unobscured cloud targets for tracking purposes (Negri and Vender llaar (1980), 
Peslen(1980), Maddox and Vonder Haar (1979)). Missing swaths of data are common with 
satellite-retrieved temperature and moisture soundings (Hillger and Vonder Haar, 1981). On 
the other hand, the Barnes scheme assumes uniform distribution of data. 

3. Tlie maximum allowable detail in the interpolated fields is to be governed by the data 
spacing. 

4. The results of the mathematical analysis of the Barnes convergence properties (appendix) 
are to be incorporated by limiting the number of passes through the data. 

5. Bounds to the grid box size should be wholly determined by the data spacing, since the 
magnitude of derivative fields like divergence and vorticity is higlily sensitive to grid size. 
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6. Measures of analysis quality should be made available to the GEMPAK user to aliow 
objective determination of the effect of user manipulation of the analysis filter control 
parameters upon the resulting objective analysis. 

7. Finally, it is essential that the GEMPAK Bames scheme be easily and quickly usable by one 
who does not possess a thorough understanding of the matliematics describing the Barnes 
scheme filter response, convergence properties, etc. appearing in this report, but wJio can 
judge the quality of the resulting objective analysis in terms of its physical content, con- 
sistency with other meteorological fields, and temporal continuity. The uniqueness of the 
GEMPAK version of the Barnes scheme rests in the manner in which these seven considera- 
tions are incorporated into an easy-to-use objective analysis package. The unique features 
are discussed below. 

a. Domain Definitions 

Unique data and mesh domain definitions were devised In order to permit easy manipulation 
of data file structures and to obtain a uniformly reliable analysis throughout the entire grid area 
displayed (Figure 3). One of these domains is termed the data file, which consists of the entire data 
set to be considered for tlie Bames objective analysis, e.g., cloud motions in the lower troposphere 
determined over a subsynoptic region for a specified time interval. 

That subset of the data file in which grid point values are computed from the data by the 
GEMPAK Bames scheme is termed the data area. Tliose observations which lie outside of the data 
area do not influence any grid point value. One conceivable data area is that in which die distribu- 
tion of the data is more uniform than tiiat within the data file. Another choice for a data area 
might be one which contains a sufficient number of observations from one data type needed to 
make a comparison with the objective analysis of data by some other type. For example, if there is 
a wish to study the impact of satellite cloud winds upon an analysis of conventional rawinsonde 
winds, then the data area should be chosen large enough to include rawinsonde data locations outside 
of the area of satellite data. 
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Figure 3. Schematie illustration of data and mesh domains used within GEMPAK Barnes scheme. 



IhQgrid display area is that portion of the entire griclded domain (i.e., data area) which is dis- 
played as “the objective analysis” to the analyst. The grid display area should lie entirely within the 
data area to avoid attempted interpretation of the analysis near the data-sparsc boundaries of the 
data area (the dashed region in Fig. 3). In this region, the data are simply extrapolated by the 
Barnes (or similar) scheme because there are not enougli data nearby to give reliable grid point esti- 
mates. A good approach to take is to have at least one or two observations lie within the data area 
existing beyond each side of the grid display area. 

b. Non-Uniform Data Distributions 

The problem of handling non-uniform distributions within the data area is presently dealt with 
in the following manner. Within any data area input to the computer, the distance between each 
observation and the nearest observation to it is determined, regardless of how uniformly distributed 
the observations are. The average of these distances over the entire data area is termed the 
“computed data spacing” An^,. Its value determines the maximum datail permitted in the objective 
analysis, However, a smoother analysis can be obtained by manually inputing any An > Ang, as 
should be done, for example, when examining the synoptic-scale impact of non-uniformly distributed, 
satellite-derived mesoscale wind, temperature, or moisture data upon vl synoptic field obtained from 
conventional data. An appropriate choice in tliis case is the An computed from the data area com- 
prised only of the conventional data. In other words, the purpose of the analysis should govern the 
choice for An, under the constraint that An > Aiij.. 

c. Control of Detail in a Two-Pass Barnes Analysis 

Once the data spacing has been defined, the analyst “fine tunes” the degree of analysis detail 
(or convergence of the interpolated field towards the observed field) by choosing a value for the 
numerical convergence parameter^ ( 7 ). In the GEMPAK version of the Barnes scheme, the maximum 


“ A sliglil degree of analysis conlrol is also made possible by manipulation of the mesh size, However, the mesh .size, 
which is allowed to vary within the limits discussed in Section 3d, should not be manipulated for the purpose of 
achieving a desired degree of analysis convergence, but rather for reasons discussed later. 


13 


detail is obtainable with 'y = 0.2, whereas the least detail results with y- J.O, for a shm value of An. 
The analyst must decide for himself how closely the final analysis should be made to fit the data. 

If one desires to see the greatest resolvable detail in the field, it is justifiable to use a small y value 
when errors in the data are only a small fraction of the signal present over the field and the observa- 
tions arc not substantially contaminated by subgrid-scale atmospheric processes. 

A limit of two passes through the data is imposed in the GEMPAK Barnes scheme (see appendix). 
Under this constraint, and the further constraint that 0.2 < 7 < 1 .0, the range of analysis responses 
permitted can be found (Figure 4). The first pass values Dq(X) are arrived at by inserting the value 
Dg(2 An) = 0,0064 for integral values of (X/2An) in (A18); the second pass values Dj^ (X) arc then 
calculated from (11), The value of 0.0064 for 0^(2 An) gives a second pass response of - e' ’ 
at the 2 An wavelength when 7 = 0.2. Since the maximum response and best fidelity cliaraeterislics 
of the Barnes low-pass filter are arrived at by choosing this 7 value, use of Dp(2An) = 0,0064 lets us 
obtain a baseline value by which responses at multiples of the 2An wave can be calculated. It is important 
to realize that under these conditions, the weight parameter is fixed by the data spacing, since 
when X = 2An is inserted into (4) with 0^(2 An) - 0.0064, 


«o = - 


'2An' 


7T 


In Dp(2An) 


= 5.052 


2An' 


( 12 ) 


In other words, the weight parameter Kq is fixed by the data spacing to give maximum response 
ofe~^ at the 2 An scale (7 = 0.2). The user of the GEMPAK Barnes scheme has the option of either 
accepting this default analysis or of making a smoother, le.ss detailed one by inputing a larger 7 value 
(0.2 <y < 1.0). Increasing 7 affects only the second interpolatioii field, which is tlic one displayed, 
SLicIi that the amplitude of the shorter wavelength features is reduced. 
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U. Objective Determination of Grid Size 

An abundant literature exists illustrating bounds on the ratio between the grid box size Ax and 
the data spacing An, such that Ax/An lies in tiie range ofabout 0.3 to 0.5 (Barnes 1964, ] 973; 
Doswell, 1977;Maddox, 1980; and others). Tlie apparent reason fortlic upper limit is tiiat since 
five grid points are required to represent vl wave, and tlie minimum resolvable v'ave is of 2 An scale, 

then Ax must he no larger than one-luilf of An (as shown in i'igiirc 5). As for the lower limit, since 


calculations of derivative (luantities like divergence ajid vorticity arc highly sensitive to grid length. 


then if such derivative fields are to represent resolvable features, one must not use a grid length that 
is much smaller than the data spacing. 


Accordingly, the GliMPAK Barnes scheme imposes the constraint that 

1/3 < Ax/An < 1/2. (13) 

The number of grid points (KX by K Y) are then determined by 


KX(KY) = 1 + (x(y) data area dimension)/Ax. (14) 

Even though the GEMPAK Barnes scheme places stringent limits on the juesh size, it remains 
versatile enough that it can accommodate round-off of the computed mesh size to convenient whole 
numbers (such as 1.0° latitude, or 10.0 km, etc.). This versatility may be necessary when making 
comparisons of objectively analyzed data to numerical model output, for instance. 


e. Quality Control Indicators 

Serious “ballooning” problems can occur when an insufficiently small number of data 
determine the value at a grid point (as discussed in section 1). Li order to avoid such problems, a 
warning flag appears to the user of the GEMPAK Barnes scheme wiienever less than 3 data values 
determine any grid point value. Such a problem is particularly likely to occur when the grid display 
area is chosen so large that it extends into data-void regions. It is desirable to caution the analyst 
tliat a better choice for the grid display area should have been made. 
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Figure 5, Determining the proper relationsliip between grid size 
and data spacing from elcmcntaiy considerations. 


The warning flag appears slightly more often than it would otherwise because of the inclusion 
of a delimiter on the computation of extremely small weights at veiy large distances between a grid 
point and a datum. Its insertion into the program increases the efficiency of its operation without 
noticeable loss of any accuracy. Presently, this cutoff distance has the value Rg = /Cg Vio, whicJi is 
4,5 times as large as the “radius of influence” R at which the weight w,„ = e“' (sec (1)). At distances 
r> Rg, the weights are assigned zero values, whereas in actuality it can be shown that the weight 
value at Rg is about 1,5 X 10~® when Ax = An/2. 

Another measure of analysis quality is the computed rms difference between the interpolated 
and observed fields. This number is displayed to the analyst after both the first and second passes, 
through the data to enable one to see the amount of rms reduction between passes and how much 
greater this reduction is with the use of a smaller 7 value. This number can then be compared to 
the known observational errors. 

Finally, computer limitations (section 4) dictate that KX and KY be no larger than 40. 


Problems arise when an analysis is attempted of clumps of closely spaced data separated by large 
distances from other clumps over a sizable data File. If one wishes to make a detailed analysis using 
An = Aiig, then several separate analyses over subsets of the data file must be made. 
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Uuspitc tlio last llinUation, tlia Gl*'Ml*AK Barnes sclicmc remains a versatile aiul easy-to-use 
objective analysis package. By incorporating a weight selection procedure based purely upon the 
data spacing, filtering the 2 An wave from the analysis, and determining the grid size from the data 
spacing, the scheme is objective in nature. Yet, by permitting user input of arbitrarily larger data 
spacing than actually exists, user manipulation of the various spatial domains, and user selection of 
the T and Ax values within bounds, the scheme retains its necessary flexibility. In the following 
section, it is shown how use of an interactive, menu-type format ami various cursor/display controls 
on a computer CR'r enhances the ease with which the GBMPAK Barnes scheme can be used. 


4. OPERATION OF THE GEMPAK BARNES SCHEME ON AN INTERACTIVE COMPUTER 

A CRT computer graphics terminal on AOIPS is utilized to permit the user to step through a 
series of menus (questions) in order to select appropriate Tiles and parameters for the GEMPAK. 
Barnes objective analysis, For each analysis the user must select 1) the vertical level for the analysis, 
2) the data sets to be included in the analysis, 3) the data and grid display areas, 4) values for the 

observation spacing An, grid spacing Ax, and 7 , and 5) the meteorological parameters to be 
analyzed. When appropriate, the values chosen most recently will be displayed as tlie default values, 

Unless default values from a previous analysis are available the following defaults are used and 
may be chosen by entering 0 on the terminal. 

(1) The default vertical level is the earth’s surface. 

(2) The default data set is the sounding data set used most recently in GEMPAK. 

(3) The data and grid display areas both correspond to the data file area. Areas are defined by 
the lower left and upper right coordinates. Since it is nccc.s.sary to align the areas to lie on 
grid points, the lower Icit corner will be fixed iiiul the upi>er rigid eoriiei moved up aiul 
right until it coincides with a grid point. AN is set to AN^; AX - AN/2. 




18 




(4) Delault values arc AN » ANj. (value cominitctl within ilcraull data area), AX - AN^./2, 
and y » 0.30. 

(5) The meteorological parameters are not defaulted; they must be selected by the user. 


When u user works In latttude/longitude coordinates, he has the capability to display the data 
and grid display areas on the CRT. Then, by manipulation of a joystick on the console, these 
analysis boxes can be easily rcsliapcd and moved until, by visual inspection, the desired result is 
obtained. After selecting which files and parameters to use, botli the actual data and the ohjective 
analysis are displayed. If the analysis is unsatisfactory to the user, he may then alter the inpu, 
parameters to the scheme by stepping through the menus once more, and subsequently see another 
analysis. Of course, instead a different parameter may be analyzed at this point. 


Following are the menus in the order in which they appear to the user and a description of the 
required responses. 

Enter Vertical Level (Default = Surface): 

A single vertical level to be analyzed, in the units of the vertical coordinate system of the sounding 
dataset name (see below), must be entered (with a decimal point included). A <CR> (carriage return) 
or 0.0 will result in default to the surface level. 

If the surface level is selected, the data set file labels will be checked to ensure that surface 
data is available in each dataset. If it is not, an error message is printed and the user may select 
another level. 

Sounding Data Set Number 1 

Enter Sounding Dataset Name (Default Disk - ) 

(0 ) or (No Dormill): 
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Up to sounding datasets^ can be entered by name. The default disk and user’s UlC (User 
Identification Code) will be added to the name if not specified 5 also “.SND” will be appended to 
the name if a file type is not entered. When a default file name appears, the file may be selected by 
entering a “ 0 ". 


An incorrect file selection will result in an error message and the file jirompt will be redisplayed. 
‘I he file labels will be checked to ensure that they have a common vertical coordinate system. 


A <CR> will end the file selection process. 

Average minimum station spacing 0 “ old values 
default == compute: 


This prompt is displayed only if default parameters from a previous objective analysis are available. 
Entering a “0” will restore the previous parameters, including the average minimum station spacing, 
AN 5 , previously computed from the present datafile, 

Enter 2 bit planes for data and grid display areas 
Default = no display: 


This prompt appears only if the user has allocated a CRT display terminal and has previously plotted 
.sounding data.*^ 


^ A GHMPAK sounding dalascl is a collection of soundings formatted through one of the options in the 
sounding input program. 

'^Each display terminal is presently a Hazeltlne terminal which has one memory of 512 X 512 X 8 bits, 
addressable as 8 bit (graphic) planes. 
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The two bit planes (1-8) must be entered separated by a comma, The first bit plane will be 
used to show the data areaj the second will display the grid analysis area, 'I’lie same bit plane n»ay 
be entered for both displays. Note that any previous displays on those bit planes will be erased. 


The menu shown in Figure 6 lists Information about the data and grid analysis areas as well as 
the current values of AN, AX, and 7 , For reference, the latitude/longitude coordinates of the 
lower left and upper right corners of the ilata file area, and the number of stations reporting data 


at the vertical level selected, are displayed at the top. Additionally, the numbc' )fsUitions within 


the data area arc displayed to tlic riglit of the station spacing, 


In 1., the selected value of AN is shown. Also listed are the last computed station spacing and 
the random station spacing. The average minimum station spacing, AN^, witliin a data area Is 
computed only if a “-I" is entered in response to this menu. The uniform station spacing is the 
average station spacing within the data area assuming that tlic stations are evenly distributed. This 


number should be used only for reference. If this number is different by a factor of 2 onnore from 
the computed minimum station spacing, then the data set is not really uniform in nature. 


2 . and 3 . show the lower loft and upper right corners of the data area, in degrees latitude, longitude. 


4. is the grid spacing, AX. Note that the grid spacing and the station (observation) spacing arc 
both in degrees latitude. The grid spacing is initially set equal to one-half of the station spacing. 

5. displays the number of grid points covering the data area in both X and Y directions. 


6 . shows the grid display coordinates in units of grid points, i.e., the subset of the data area 
I'.iitl wliidi will beanm’ llic gritl disiday area. The laliliide. loiigidule comers of llie j'i'ul dispt.iy aivii, 
which cannot be altered directly, are shown under the coordinates. 
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♦^♦BARNES analysis MENU*** 



(28.8 -106.0 

41.4 

-86.6) 

34 STATIONS 

1. 

STATION SPACING 


1.93 

32 STATIONS 


LAST COMPUTED SPACING 

1,93 




UNIFORM STATION SPACING 

2.E1 




-1 TO COMPUTE STATION SPACING 



2. 

DATA AREA (LOWER LEFT) 

28.82 

- 

-106.37 

3. 

(UPPER RIGHT) 

41.36 


“86.11 

4, 

GRID SPACING 

0.965 



5. 

# GRID POINTS (KX, KY) 

21 


14 

6. 

GRID DISPLAY COORDINATES 

1,21 


1,14 


(LOWER LEFT) 

28.82 

- 

106.37 


(UPPER RIGHT) 

41.36 


-86.11 


Figure 6. Example of GEMPAK Barnes “menu” of input parameters. 





Hnter number and new values 
0 = initialize P = PRINT D = Data Box 
Def£.ult - accept: 


This prompt appears following tlie menu shown in Figure 6. The user may enter a menu selection 
number (1-6) followed by a comma and new parameter value(s) separated by commas, The values 
entered are checked for obvious errors, changes are made to the required parameters and the menu 
is redisplayed. A station spacing AN > AN^. can be entered at this step to account for non-unifonn 
data distributions. Entering a new value for one parameter can affect values for the other parameters 
(e.g., new station spacing alters data area slightly). 

A response of “0” will reset the menu to the values initially displayed. “P” will print the menu 
on the line printer. “D” will display a yellow box on the display screen. The box may be moved 
and shaped to redefine the data area using the “MOVE” and “SHAPE” buttons on the terminal 
console. The “DEFINE” button is used to accept position of the box. The user may enter “:XF” 
to turn the box off. 


ENTER GAMMA BETWEEN 0.2 and 1.0 
DEFAULT GAMMA = 0.3: 


Any value for 7 within the limits 0.2 < 7 < l.O is acceptable. If any other value is entered, the 
menu is displayed again. A <CR> will default to 7 0.3. 

Enter list of parameters to be analyzed (< = N): 
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Up to N parameters (e.g., temperature, wind speed, etc.) may be selected by number from the list. 
N is limited by computer considerations. Since the most time-consuming part of the objective 
analysis is the exponential computation, large grids or data areas containing a large number of 
stations will limit the number of parameters that can be analyzed at one time. The number of 
exponentials computed is approximately 2 X M X KX X KY, where M is the number of stations. 
Since thi.s numbci is not a function of the number of parameter.s being analyzed, It is desirable to 
compute as many parameters as possible during a single analysis. 


PASS 1 

PARAMETER 2 

RMS = 2,38 

#STATIONS = 32 

PASS 2 

PARAMETER 2 

RMS = 0,74 

# STATIONS = 32 


This message prints the RMS after each pass. No user response is required. The reduction in RMS 
error between the two passes is greater when a small 7 value is used. The RMS at the second pass 
can be compared to known observational uncerhiinties to help the user make a judgment of analysis 
quality. 


Following each pass, parameter values are computed at each station using bilinear interpolation 
with the four surrounding grid points. The RMS is computed using the differences between tlie interiio- 
lated station value and the actual station value (see ( 8 ) in Section 2). The parameter number refers 
to the GEMPAK number for the parameter being analyzed. The number of stations actually used 
in the analysis is also shown. 

THERE ARE POINTS IN THE GRID DISPLAY AREA WITH 

INSUFFICIENT DATA. 

A LISTING WILL BE SENT DIRECTLY TO THE LINE PRINTER UNLESS 
AN “N” IS ENTERED: 
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The prompt is displayed only if fewer than three data values have influenced any grid point, as 
mentioned in Section 3e. 


ENTER GRIDDED DATASET NAME (DEFAULT DISK = ) 
(0 = ) OR (NO DEFAULT) 


Tliis gridded dataset, a file which stores up to 20 grid analysis arrays, and has the form “NAME GRD,” 
must be entered. Again, the default disk and the users UIC will be added if not specified here; “GRD” 
will be appended if a file type is not entered. When a default file name is given, the file may be 
selected by entering a “ 0 ”. 

OUTPUT TEMP IN GRID NUMBER 10 
ENTER NAME FOR GRID: 


This prompt asks the user to specify a name for grid 10; in this case, the gridded data resulting from 
an objective analysis of temperature (TEMP) is stored in grid 10. A 12-character name may be en- 
tered to identify the grid in the file. The grid number is the location of the grid in the file. The 
grid is accessible for purposes of algebraic manipulation (e.g., pressure tendency calculations made 
by subtracting two gridded pressure datasets) and contouring. 

5. METEOROLOGICAL APPLICATIONS OF THE GEMPAK BARNES SCHEME 

In this section, station and grid spacings (An and Ax), data and grid display areas, and the 
numerical convergence parameter ( 7 ) are each varied to examine the resulting effect upon the 
objective analysis of subsynoptic-scale wind data. Each of the resulting GEMPAK Barnes analyses, 
except for the last one to be discussed, will be of a uniformly distributed wind data set obtained from 
rawinsonde data collected during the SESAME (Severe Environmental Storms and Mesoscale Experi- 
ment). The last analysis will be that of non-uniformly distributed wind data composed of both the rawin- 
sonde observations and satellite-derived cloud motion estimates. The purpose of the latter analysis is to 
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show how An can be chosen l.’rger than the computed value An^ so as to obtain an analysis more 
representative of the scale resolvable by the uniformly distributed, conventional data set. 


The data file used for all but the last analysis is derived from the SESAME regional-scale 
network, a uniformly distributed array of rawinsonde observations including 23 National Weather 
Service stations and 19 siipplcmental sites over the .south central United Stales (h'igure 7), In 
particular, an analysis of wind speed at the 310 K isentropic surface at 2100 GMT on April 10, 1979 
is made from data existing at 32 of the 42 possible stations. These 32 stations constitute the data 
file^ wliich form the base for the objective analysis. The latitude and longitude bounds of the data 
file are: 


Latitude Bounds 
upper: 41.36 
lower: 28.82 


Longitude Bounds 
western: -106.37 
eastern: -86.11 


(15) 


a. Map Generation 

The first procedure in generating a Barnes analysis is to display a map of the obsemtions 
within the data area (the only real requirement is for the analyst to display just a map, but the 
inclusion of the observations could minimize errors related to the placement of data and display 
areas). A map containing the SESAME wind speed observations (Figure 8) is generated and displayed 
by GEMPAK’s sounding plotting program. For this example, the coordinate bounds of the map were 
chosen as follows: 

Latitude Bounds Longitude Boujids 

upper: 43.4 western: -108.4 (16) 

lower: 25.4 eastern: -83.4 

These coordinate choices produce a map covering 25° of longitude and 18° of latitude, whici) incliule.s 
all the observations within the data file plus an additional surrounding data-free area of about 3°. 
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As long as no furtlier changes are made to the map’s coordinates within the sounding plot program, 
the latitude/longitude coordinates of the map area will be retained throughout tJie use of the Barnes 
program, so that the data and display areas can be accurately placed over the map. 

b. Station/Grid Spacings and Data/Grid Display Areas 

Once tlie map has been generated, the GEMPAK Barnes menus are utilized. Following the 
procedure outlined in Section 4, the user first specifies liis desire to analyze the wind speed at verti- 
cal level = 3 10 (in this case, 3 1 0°K), He is afterwards faced with the task of generating suitable data 
and grid display areas. 

The interactive CRT tenninal allows the analyst to create and alter two “boxes” on tlie CRT 
display screen, one box enclosing the data area, and the other enclosing the grid display area. During 
the initial internal computation of Ax and An, only one box will appear since the program will 
assign the same coordinates to both the data and grid display area. Unless old values for An arc 
available, the average minimum station spacing v/ithin the data file will be computed, which in the 
present case (Figure 8) is An^ = 1.93°. Under the constraint (13), the grid spacing is initially set to 
its maximum allowable value of Ax = An^/2 = 0.965°. Upon input of the data file bounds (15), 
the values of KX and KY are computed internally from (14), and hence the Barnes Analysis Menu 
will appear on the CRT display screen with the parameter values (shown earlier) in Figure 6. 

The calculated value of An^ may be modified now to meet the needs of the analyst. It is essen- 
tial to realize tlie importance that such a modification has upon the analysis (see Section 3c), within 
the ’’mit that any input An must be no smaller than An^. For discussion purposes, the value is in- 
creased slightly to An - 2.00° in order to allo w an internally calculated grid spacing of Ax = 1.00°, 

In general, the value of An should not be arbitrarily increased beyond a few percent of its calculated 
value unless a non-uniform data area is to be analyzed, for otherwise details resolvable by the data 
will be lost. 
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TJje change of An from 1,93° to 2.00° results in a slight enlargement of the box representing 
the data area, While holding the lower left corner of the original coordinates (15) fixed and 
accounting for the grid spacing Ax of 1.0 degree, the GEMPAK Barnes program computed new co- 
ordinates for the data area: 

Latitude Bounds Longitude Bounds 

upper: 41.82 left: -106.37 (17) 

lower: 28.82 right: -85.37 

Note that the coordinates of the upper right corner of the data area have changed to allow the box 
to increase in size while including all observations from the data area, Tiie degree differences 
between the latitude bounds and the longitude bounds are now multiples of the grid spacing Ax of 
1 .0 degree. 

Following the clioice of Ax, the data area is chosen by either typing in the coordinates or by 
using a joystick-controlled function tliat shapes a rectangular box around the data area. It is easier 
to position a box manually than to determine the coordinates of the lower left or upper right edges 
of the data area. In either case, once the data area has been entered, the program may alter the co- 
ordinates and the box slightly to accommodate the grid spacing. Tiie resulting data area box is the 
larger of the two boxes in Figure 9. 

The grid display area should be entered as a subset of the data area. A good approach to take 
(see Section 3a) is to allow at least one or two observations to lie beyond eacli side of the grid 
display area. Since the slight increase of An from 1 .93° to 2.00° results in an increase of the grid 
array from 21X14 (Figure 6) to 22 X 14 when Ax = 1.00° (see (17)), a grid display area 
imbedded two grid points (or 2°) within each side of the data area requires the following sequence 
of valtK's fo be inpul in item #6 of Ihe names menu; 

3 20 3 12, 
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Figure 9. Data and grid display areas selected for analysis of SESAME data set (see Figure 8). The 
boxes were obtained both by typing in the coordinates and by using joystick-controlled 
functions on the AOIPS computer console. Note size of grid mesh in lower left corner. 
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iiml rt'sulls in an 18X10 array of grill poini.s williin a jt»m/ ilisplav ami with the roltowing 
coordinates: 


Latitude Bounds Longitude Bounds 

upper; 39.82 west: *-104.37 

lower: 30.82 east; -87.37 

The grid display area is shown enclosed by the smaller of the two boxes in Figure 9. 

A comparison of two objective analyses of the SESAME data set generated by (1) having a grid 
display area smaller than the data area (Figure 10a) and by (2) having the two areas represented by 
one box (Figure 10b) demonstrates the rationale behind choosing the former option. A less reliable 
analysis results when the data and grid display areas coincide (Figure 10b), especially near the cast 
and south sides of the box, Notice that the observations at Nashville, Tennessee and througliout 
southern Texas are not incorporated into that analysis, resulting in loss of valuable information on 
the strong wind speed gradients actually there. On the other hand, the analysis displayed in Figure 
10 a, which was generated from values interpolated to grid points enclosed by the outer box (the data 
area) but displayed only within the inner box (the grid display area) is more reliable because observii- 
tions exist in the region between the grid display area and the data area. Therefore, each grid point 
within the grid display area has been influenced by values at station locations which totally mrroimd 
it. 

c. Numerical Convergence Parameter 

In the case of the two analyses just described, the numerical convergence parameter 7 was 
assigned a value of 0.3. The GEMPAK Barnes program allows 7 to range in value between 0.2 and 
1.0. Analyses of tlie SESAME wind speed data set with four different values of 7 (0.2, 0.3, 0.5, and 
1 .0) are displayed in Figure 1 1 to Illustrate the impact of 7 variations Upon analysis detail (recall 
Section 3c), Tiie analysis produced with 7 = 0,2 exhibits significant detail at the 2 A n^ scale. As 
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GAMMA « .6 



Figure 1 1. CJianges in the objective analysis of the SESAME wind speed data set 
(see Figure 1 0) brought about by variations in the value input for 7 . 
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discussed in Section B of the Appendix, tins detail is likely to he the undesirable result of aliasing and 
of random erroi^ in the data, Larger values of 7 lov/er the response to that scale (as given by ( I • ) and 
Figure 2) and accordingly dampen the analysis detail. 

Tlie analyst must judge which value of 7 produces the best analysis for his purposes. Altliough 
lliere is no intenlion in Ihis rept>rt lo coverall aspects of (he general I'aciors iiivolveil in ohiective 
analysis quality, a few points should be kept in mind. First, the nm iliffercntu between the inter- 
polated and observed fields should be much less than the estimated rms error in the observations; 
tliis guideline should help serve as an upper limit on the choice of 7 , Second, the meteorological 
features appearing in the analyses should exhibit acceptable temporal and spatial continuity. The 
lower limit on 7 is determined by the continuity considerations, which is best made by examining 
derivative fields like divergence and vorticity. Third, the analyst should examine tlte interpolated 
fields resulting from the largest 7 value first, and then reduce 7 in steps until an unacceptably small 
signal-to-noisc ratio becomes evident. 

Application of the rms quality criterion to the objectively analyzed wind speed fields in Figure 
11 is illustrated first. An observational rms error of 3 to 6 m s“' can be expected for winds between 
the 700 and 400 mb pressure surfaces, through which the 3 10 K isen trope passes in this case {Fuel- 
berg, 1974), The computed rms differences between observed and objectively analyzed fields ex- 
ceeded 1 .5 m s”^ for 7 values equal to or greater than 0.5. Thus, tlie 7 = 0,2 or 7 = 0.3 analysis 
would be selected upon this basis. 

The other measures of analysis quality can be derived from the calculated velocity divergence 
fields at 3-hourly intemls from 1500 GMT to 2100 GMT on this date. The fields shown in 
Figures 12, 13, and 14 are generated from the objectively analyzed u and v wind component fields 
using 7 values of 0.2, 0,3, and 1.0, respectively. Examination of thc.se fields for spatial and 
temporal continuity is made, bearing in mind that, althougli the archived SESAME data have been 
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Figure 13. Same as in Figure 12, except 7 = 0.3. Labelled features were followed from one analysis 
to the next. 
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Figure 14. Same as in Figure 12, except 7 = 1 .0. 











checked for gross errors, stricter quality control measures were not applied to the data set to 
account for any possible bias errors. 

Each set of analyses exhibits at least some continuity, with the degree of continuity decreasing 
as y is reduced. Features in the 7 = 1 .0 analysis are easy to follow, but lack detail, which is consistent 
with its large rms error value. Many more features appear in the 7 = 0.2 analysis, but the fact that 
their amplitudes and propagation velocities fluctuate greatly indicates an unacceptably low signal-to- 
noise ratio in the analysis. This observation is clearly evident in the case of the isolated convergence 
maxima (-3 X 10”^- s"^ ) over eastern Oklahoma and western Missouri at 1800 GMT (Figure 1 2b). 

Such small-scale features are evidently not resolvable with this particular data set. 

In contrast, the alphabetically labelled features in the 7 - 0.3 analysis can all be traced with 
acceptable temporal continuity over at least two time periods. Furthermore, these features move at 
velocities (20-30 m s“^) characteristic of the observed wind speeds (see Figure 8) and in a northeasterly 
direction consistent with the prevailing advecting flow. Other unlabelled centers either dissipate 
between periods (as with the divergence maximum over southwestern Arkansas at 1 500 GMT) or move 
out of the gridded display area too quickly to be followed (as witln the divergence maximum over the 
Texas Panhandle at 1500 GMT). Thus, the 7 = 0.3 analysis provides the best results in terms of rms 
error and temporal/spatial continuity criteria. 

d. FTfect of Grid Size on the Objective Analysis 

Recalling that for a given data spacing An, the GEMPAK Barnes program constrains Ax to 
range between An/3 and An/2, then for the SESAME wind example whereby An = 2.0°, Ax must 

lie between 0.65° and 1.0°. The latter case appeared in Figure 10a, and for comparison the case 
Ax = 0.65° appears in Figure 15. The differences between the two analyses are insignificant. This 
comparison cmpliasizcs tiic fact that vanulioiis of llw grid .size wiiliiii Hie ullowahie iiinii.s ha.s u 
negligible effect upon detail in the objective analysis, and that it should not be varied with that 
purpose in mind. Rather, it should be varied for such purposes as to ease the discussion (as here, 
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Figure 15, Effect upon objective analysis of reducing grid size to its smallest value permifed within 
the GEMPAK Barnes scheme (Ax = An/ 3), This analysis should be compared to that in 
the upper portion of Figure 10. 
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An was increased from 1 .93® to 2.00° to make Ax = 1 .00° a convenient round number), or to make a 
one-to-one comparison between objectively analyzed data and numerical model data, 

e. Objective Analysis of Non-Uniformly Distributed Data 

A comparison is made between an objective analysis of SESAME rawinsonde reported winds at 
825 mb at 1800 GMT on April 10, 1979 (Figure 16a) and an objective analysis of those same winds 
when combined with non-uniformly distributed satellite data (Figure 1 6b). These analyses are being 
used in an on-going investigation at NASA/GLAS to determine to what extent ai signment of the cloud 
motion vectors to an incorrect height will degrade a conventional analysis of winds at the SESAME 
regional scale (Peslen, et ah, 1982). With this purpose in mind, the grid display area is adjusted to 
just barely cover the area of satellite data. Likewise, for this purpose, the value of An input to the 
computer is that of the uniformly distributed, conventional rawinsonde data of the SESAME 
regional scale (An = 2.0°), rather than the computed minimum observation spacing A n^, which has 
a very small value because of the small distances between adjacent satellite observations. Moreover, 
use of An^. would have been impossible on the AOIPS computer (see Section 3e), because Ax would 
also have been too small and hence KX (KY) too large. 

Clearly, the satellite data do exert an impact upon the conventional wind analysis on the scale 
resolvable by the SESAME regional rawinsonde network. The question remains, however, whether 
the satellite data can provide useful meteorological information on scales smaller than is resolvable 

by these rawinsondes. At least part of the difference between Figures 16a and 16b can be ex- 
plained by a systematic mis-assignment of all the cloud vectors to one isobaric level, when in 

actuality the clouds exist at different heights above the local terrain (Peslen, et al., 1982). The 
problem of mis-assignment of heights is most crucial in a vertically sheared environment. The 
following questions need to be unravelled; (1) what coordinate system to assign cloud vectors to 
with minimum error; (2) what the smallest scale is at which satellite wind data provides reliable in- 
formation; and (3) what observation spacing An is the proper one to employ with such non-uniformly 
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Figure 16. Comparison of objective analyses of uniformly (upper figure) and non-uniformly (lower 
figure) distiibuted wind data sets. Data spacing in lower analysis is set equal to that in 
the upper analysis in order to examine impact of satellite-derived cloud motion data 
upon scales resolvable by the rawinsondes. 
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distributed data. Further refinement of such anu^yses as tliat in Figure 16b can only be made follow- 
ing the unraveling of these questions. 

6. SUMMARY 

The specific adaptation of tlie Barnes (1973) objective analysis scheme to GEMPAK on the 
AOIPS computer at NASA/Goddard has been described in this report. The Barnes scheme was 
selected because it is computationally efficient, has filter response characteristics that are known 
functions only of the data, and adequately recovers details on the second pass through tlie data 
even when a large influence radius is used to insure that sufficient data influence is exerted at all grid 
points. A mathematical analysis proves that the scheme is absolutely convergent, and that only two 
passes through the data are needed to enable the analyst to control resolvable small-scale details in 
the interpolated fields. 

The GEMPAK version of the Barnes scheme is unique in that it makes possible a sound 
objective analysis of any two-dimensional scalar data field by incorporating several objective con- 
straints upon the analysis over whicli the user-analyst has no control, while still providing .some ii.scr 
input via an interactive computer video terminal. Thus, the GEMPAK Barnes scheme is altogether 
versatile, objective, and practical (Figure 17). 

By permitting user selection and alteration of the following input parameters, the scheme re- 
tains its necessary versatility : 

(1) Data area as subset of a data file. 

(2) Gridded display area as subset of the data area. 

(3) Data spacing An. 

(4) Grid spacing Ax. 
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Figure 17. A summaiy description of tlie unique features of the GEMPAK Barnes sdieine 


(5) TJie numerical convergence parameter 7 , within tiie limits 0.2 < 7 < 1.0 imposed 
(determines the degree to which tlie interpolated field converges to the observed data 
field, i.e., the amount of analysis detail). 

By retaining the following objective features, which are unique features of the GEMPAK Barnes 
program, a reliable objective analysis can be made without unnecessary subjective human intervention; 

(1) Tlie detail in the analysis is constrained by the data distribution as the only permitted An 
fliat can he inscrled is t)iiu larger than An^. (computed minimum spacing l)ctwecn data 
points, averaged over the data area). 

(2) Bounds are placed upon Ax by the data spacing, namely An/3 < Ax < An/2, to insure 
proper represontatign of resolvable features. 

(3) The weights of the Barnes low-pass filter are determined solely by the data spacing, and 
are calculated internally with no user control once An has been input. The most detailed 





analysis permitted is that one resulting from input of An = An,, and y « 0.2, in which 
case the 2 An^ wave is assigned a mere 1 /e response amplitude. 

(4) No user control of the number of passes througli the data is permitted, as this number is 
fixed at two based upon a mathematical analysis of the convergence properties of the 
Barnes scheme (presented in the Appendix). 

Use of an interactive, menu-type format and various cursor/ display controls on a computer 
video screen makes the GEMPAK Barnes scheme an easy and practical operation. The user ^^'’ps 
througli a series of menus (questions) displayed on the screen to help in the selection of input 
parameters, quickly sees the data and grid display areas and the contoured, interpolated fields pre- 
sented on a map, and can then alter the values of the parameters to test the effects upon the fields. 

In addition, several quality control indicators are presented to him ('such as rms difference between 
the interpolated and observed fields). Practicality Is augmented by being able to use the same calcu- 
lated weights for many parameter fields. Thus, a person with only a superficial knowledge of the 
Barnes scheme can quickly obtain sound objective analyses that faithfully represent the data. 

Two cases were presented to illustrate the versatility, objectivity, and practicality of the GEM- 
PAK Barnes scheme. The first, that of a unifonnly distributed rawinsonde data set, is one in which 
An was assigned its calculated value An^, so that the 2 An^ wave was essentially filtered from the 
gridded fields. It was demonstrated that when 0.2 is selected for the y value, the maximum detail 
for a given data spacing is realized in the gridded fields. Such a small y value is justifiable only when 
the data are not substantially contaminated by errors and subgrid-scale atmospheric processes (which 
may unrealistically alias energy to larger wavelengths). In the case presented, the data were rather 
good and y = 0,3 was chosen as the best analysis. The effects of varying the other input parameters 
upon the gridded analysis were also shown. 

In the other case, inclusion of satellite-derived cloud motion vectors with the rawinsonde data 
resulted in a non-uniform data set. Use of An„ would be inappropriate here, because it has yet to 
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be determined what the smallest scale is at which such satellite data provide reliable infomiiation. It 
is for such a reason as this that the user is allowed to intervene (by inputting a larger An than the 
calculated one *, with the result that a less detailed analysis is made. 

The practicality of the GEMPAK Barnes scheme is demonstrated by the fact that each of the 
analyzed maps in the cases just presented were produced in hard copy form at intervals of 3-4 min- 
utes once the set of input parameters had been decided upon. For an experienced AOIPS computer 
user, it takes one 1-2 hour session to generate the map display and go through the variations leading 
to the final decision of which parameter values to use in future map analysis generation with the given 
damset. 


'7, FUTURE PLANS 

Several modifications to the GEMPAK Barnes scheme are botli anticipated and possible with 
present AOIPS computer resources. Upstream-downstream enhancement of the weighting parameter 
•Cq would be preferable in jet-like regimes where isolines of the analyzed scalar tend to align with the 
v'Ind direction (Endlich and Mancuso, 1968). Use of the following modification 

= KqH+^cos^ 4>) (19) 

has been suggested by Barnes (1973), where |3 -V/V^ <p is the angle between the v/ind vector V and 
the vector that points between grid point and observation, and V*^is a scaled wind speed with value 
ot from 1 to about 3. In the case of hurricane observations, a modification of to account for its 
circular flow (as attempted by Bergman and Carlson (1975)) is desirable, perhaps based upon the 
Ranldne vortex model. 

One highly desirable addition to the present package is a comprehensive, automated, interactive 
data editing/quality control routine that would allow tlie analyst the power to selectively alter or 
omit troublesome observations. Displays of data time series, time tendencies, map plots, hydrostatic 
and superadiabatic checks, etc, are envisioned for this addition. 


Additional work is needed to more objectively determine what data spacing to use with non- 
uniformly distributed data sets. It is likely that the optimum An will be a function of the repre- 
sentativeness, or error magnitude, characteristic of tlie particular data type Jis well as the nature or 
degree of non-uniformity of the data distribution. Map displays of botli the influence radius and flic 
difference between observed and interpolated values at actual station locations would help one sec 
the spatial variation in accuracy of the analysis due to non-uniformity of the data distribution. 

Access to a larger computer can make possible several other extensions to the core GRMPAK 
Barnes scheme: (I) use of time series data can enhance detail in the analysis by adding more data 
in space through a time-to-space conversion process (Barnes, 1973); (2) a three-dimensional Barnes 
analysis scheme would better serve the needs of dual Doppler radar meteorologists than a two- 
dimensional one;^ (3) finally, it is envisioned that two-dimensional gridded fields produced by tiie 
GEMPAK Barnes scheme can be used as initial fields for three-dimensional mesoscale analysis 
schemes that require “data at grid points,” such as variational schemes that incorporate various pliysi- 
cal constraints into tlie analysis. 
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APPENDIX 


The mathematical-numerical analysis presented here shows the Barnes (1973) scheme is abso- 
lutely convergent, but that no significant improvement in filter response fidelity (steepness of 
response curve) is obtained by making more than two passes through the data when 7 is chosen 
sufficiently small. It is concluded that only one correction pass is sufficient to achieve acceptable 
analysis convergence at resolvable wavelengths, 

A. Convergence Proof 

Convergence in an absolute sense occurs when the difference between the observed and analyzed 
data fields vanishes as the number of passes through the data approaches infinity, or stated mathe- 
matically 

lim [f(x,y) -gN(x,y)i ~0, • (Al) 

N-^oo 

where gj,^(x,y) is the interpolated field obtained after N iterations (N+1 passes). It was shown in 
the main text that a single application of the reduced weight parameter /Cj (2) upon the weight function 
Wj^(l) results in the true response function D*(l 1) and the “difference field response function” 

Dj (6) at the second (“correction”) pass. Weight functions at additional passes using successively 
decreasing values of assuming that 7 is kept constant through all passes, are defined by: 

Xn = 7Kn-1 "" (A2) 

Application of these additional filter functions resuiv “difference response functions” at each 
pass given by 

Dn = , <A3) 

When the value chosen for 7 is less than unity, the response is further accentuated at each additional 
pass, particularly at short wavelengths where the initial response is small. 

It follows from (5) that the third pass (N=2) interpolated field is: 
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g2=gi+(f-gj)D2, 

or upon substitution of the expression for gj from (5): 

§2 = 80 + (‘*■"80) [Di + D2(l-Dl)]. 

Following tliis procedure one step further, one can easily shew that on the fourth pass: 

83 = 82 + (f-82)t>3 

' = 80 + (f - 80) IDi + D2(1 - Dj) + Dgd - D2) (1 - Dj)] . 

Accordingly, the general form of the equation describing the Nth pass interpolated field is: 


(A4) 


(A5) 


(A6) 


f N 1 V 

:^ + (f-g„) S [Dj n (1-D. 1 )] . 

\ i=2 j==2 ) 


= + + 2 ; m 11 U-U._dJJ. (A7) 

To obtain the condition for convergence, this equation is substituted into (Al), resulting in: 


( N i 

lim (f-g^){l- D] + 2 [Dj n (1-D: = 

N->~ I I i-2 j=2 }) 

Finally, because (f - go) is a fixed constant, convergence is attained when 


(A8) 


'/'N “ ^ j 


(A9) 


where is given by (6) and 


N i 

2 [Di n (l-Dj,,)] 
i=2 j=2 


2<N<~ 


(AlO) 


Before proceeding with the analysis of the convergence criterion, it *s important to understand the 
relationship between it and the actual response function at the Nth iterative pass D^, defined by 

liiHi > 


D' 


N “ f • 


^.All 
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Use of (A7) in (A1 1) results in the simple relationship: 


D5 = Do + (l-Do)(Di + i//f^). 2<N<oo. (A12) 

Thus, convergence defined in the equivalent sense to (Al), namely 

lim gM = f , (A13) 

demands that D^, as defined by (A12), approach unity in light of criterion (A9). 

The nature of the convergence properties of the Barnes (1973) objective analysis scheme can be 
understood from an analysis of (A9), (AlO), and (Al2). It can be shown analytically that (AIO) is 
a convergent power series, and numerically that it converges to the value given in (A9). Applying 
the ratio test for convergence to (AIO), we have 

N+1 

Dn+ 1 n (i“Dj_p 

i=2 ^N+1 

i = (A14) 

Dn n (I-Dj.i) 

i=2 

The series (AIO) is absolutely convergently only if L<1. Substitution of (A3) into (A14) results in: 

L = (A15) 

Recalling that is constant for any given choice of weight factor and has a value 0 <Dq < 1, 
then L must be a constant, as required for convergence. In the limiting case as N-^-oo, both 
and 7^->0 for the range 0<7< 1 , so that L->-0. Thus, there is absolute convergence in the limiting 
sense as 

The results of calculating Dj^(Dq, 7,N) on a programmable desk computer (Figure Al) indicate 
that it converges rapidly to unity, as required. Barnes (1973) noted that the fastest restoration of 
small wavelength amplitude suppressed in the first pass filtering-interpolation (D^) is obtained with 
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the smallest values of 7 (Figure 2). For a choice of 7 = 0 . 20 , the value of approaches 1 .0 to 
within ten decimal places by the N“ 6 th iteration, whereas the same degree of convergence is 
reached by N=9 when 7 = 0.45 (Dp = 0.0064). Convergence is obtained relatively quickly even for 
7 = 0.8 and is attainable in a finite time as long as 7 < 1 . Thus, the Barnes (1973) objective analysis 
scheme converges to the value specified in (A9) when 0 < 7 < 1 . It is noted that tliis scheme forces 
the interpolated field to converge much more rapidly to the observed field tlian does the Barnes 
( 1 964) schcinc, as can be easily seen by comparing (A 1 2) to equation (20) of that paper. 

B. Effect of Multiple Passes Upon Filter Fidelity 

Having thus verified that convergence is attainable, the question of whether making more than 
two passes through the data can effect a significant enhancement of the small, but resolvable, waves 
must be answered. For most purposes, it is desirable to suppress the response to a wave whose scale 
does not exceed twice the average minimum data spacing (X<2An). Considering the Gaussian nature 
to the response function at the first pass given by (4) and illustrated at additional passes by Figures 
2 and Al, then the final response at this minimum resolvable scale should be limited by 
D^(\ = 2An) «^e"^ . Under this constraint, high frequency “noise” generated by random errors and 
energy aliased from shorter wavelengtlis to larger wavelengths will be effectively filtered from the 
analysis. The aliased energy can result from both the data discretization process and the sampling 
of the atmosphere when it exhibits such subgrid-scale events as thunderstorms. 

Thus, the question can be replirased as, for a given fixed weight parameter k^, whether the 
steepness of the filter response curve at wavelengths larger than 2An can be appreciably increased 
(filter fidelity enhanced) by making more than two passes under the constraint that the final 
response at X=2An does not exceed e“^ . A comparison was made between the response curves 
generated by a two pass, small 7 filter and those generated by multiple pass, larger 7 filters, since 
the best filter fidelity is obtained when 7 is smallest (Figure 2), 

The results of the comparison (Figure A2) were obtained by employing the following procedure. 
First, the 2An wavelength was defined in terms of known quantities. Since the weight parameter 
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A/2An 


Figure A2. Effective enhancement of filter response produced by making multiple passes, as 
function of resolvable wavelengths X>2An. Each response curve results from holding 
7 constant for N iterative passes under the constraint that D^(2An) = e“^ , from (A 12), 
where X = 2An is defined by (A16) with Dg(\ = 2An) = 0.0064, and D^(X=?^2An) is 
calculated from (A18), 
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Kq must be held constant for the purpose of the comparison, and under tliis condition a singular 
relationship exists between any wavelength X and the first pass response at that wavelength Dg(X), 
then a useful definition according to (4) would be: 


2An ” n-yj-Kf^ I In DQ(2An) . 


(A 16) 


The actual value of Kq need not be known to use tliis definition for 2An, since when is constant 


JL / lnP^(2An) 
2An “V InDo(X) 


(A17) 


Thus, the effectiveness of making additional passes at larger 7 is examined as a function of multiples 
of the 2 An wavelength. Each curve in Figure A2 is the result of the search for that value of 7 which 
gives DjI^(X=2An)=e~^ (N>1), and then finding Dj^ at other wavelengtlis by the following formula 
derived from (A1 6), 

D^(X) = [DQ(2An)](2^"A)^ (A18) 


for integral values of (X/2An), inserting the result in (A12) for the chosen 7 value, and plotting 
against multiples of the 2An wave according to (A17). The selected value of DQ(2An)=0.0064 used 
in (A18) is that one wliich gives a second pass response of D^(2An)=0,37=e”^ when 7=0.2 (see 
:‘igure Al). The value of 7 =0,2 is chosen to represent the two pass, small 7 case. 

The results in Figure A2 show that when N=I,2, and 3 iterations through the data are made, 

Jie responses are 0.84, 0.92, and 0.96, respectively, at twice the minimum resolvable scale. Differ- 
ences between these responses are no more than 12% and decrease at larger wavelengths. These 
differences are further attenuated when an even smaller 7 value case is chosen as the basis for com- 
parison with multiple pass filters. Thus, one correction pass using a small 7 value provides for a 
highly acceptable degree of filter fidelity. Therefore, if one wishes to make the final objective 
analysis fit the data as exactly as possible, the same result can be obtained by making a greater 
number of interpolation passes (which is computationally wasteful and can cause greater “ballooning” 

effects in data-sparse areas), or by using a small 7 value to reduce the coiTection pass influence radius, 
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